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ABSTRACT 


We  present  a  theoretical  framework  that  can  be  used  to  treat 
approximation  techniques  for  very  general  classes  of  parameter  estimation 
problems  involving  distributed  systems  that  are  either  first  or  second  order  in 
time.  Using  the  approach  developed,  one  can  obtain  both  convergence  and 
stability  (continuous  dependence  of  parameter  estimates  with  respect  to  the 
observations)  under  very  weak  regularity  and  compactness  assumptions  on  the 
set  of  admissible  parameters.  This  unified  theory  can  be  used  for  many 
problems  found  in  the  recent  literature  and  in  many  cases  offers  significant 
improvements  to  existing  results. 
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1.  Introduction 


An  important  class  of  scientific  problems  related  to  the  control  of  flexible 
space  structures  entails  the  estimation  of  parameters  in  distributed  or  partial 
differential  equation  models.  These  inverse  or  parameter  identification 
problems  arise  in  several  contexts.  A  rather  obvious  class  of  such  problems 
involves  the  development  and  identification  of  physical  models  for  flexible 
structures  -  e.g.,  see  [BC1],  [BC2],  [BCR],  [BR2J,  [BPR],  [BWIC],  [BFW],  In  these 
problems  one  typically  investigates  models  from  structural  mechanics  which 
contain  parameters  representing  elastic  properties  (such  as  stiffness,  damping)  of 
materials.  The  inverse  problems  then  consist  of  estimating  these  parameters 
using  data  obtained  from  observations  of  the  system  response  to  dynamic  loading 
or  displacement.  The  partial  differential  equation  models  usually  are  second 
order  in  time,  higher  order  in  space  (for  example,  the  Euler-Bernoulli  or 
Timoshenko  beam  models  and  their  higher  dimensional  analogues  for  plates,  etc.). 

The  related  inverse  problems  can  prove  quite  challenging,  but  arc  extremely 
important  as  a  precursor  to  and  integral  part  of  the  design  of  control  strategies 
for  space  structures.  • 

\ 

A  second,  less  obvious,  class  of  inverse  problems  that  are  important  in  the 
study  of  space  structures  involves  those  arising  in  material  testing,  in  particular, 

in  the  area  of  nondestructive  evaluation  of  materials.  These  problems,  such  as  _ 

those  related  to  the  detection  of  structural  flaws  (broken  fibers,  delaminations,  ~~~" 

cracks)  in  fiber  reinforced  composite  materials,  have  been  addressed  at  NASA 

(for  example,  in  connection  with  evaluation  of  the  solid  rocket  motors  for  the  — 
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Space  Shuttle)  using  a  variety  of  techniques  (acoustic,  thermal,  etc.).  One 
technique  [HWW],  [BKO]  entails  the  use  of  thermal  diffusivity  properties  (a 
thermal  tomography)  to  characterize  materials.  The  resulting  inverse  problems 
involve  using  boundary  observations  for  the  estimation  of  thermal  diffusion 
coefficients  and  interior  boundaries  of  the  domain  (i.e.  the  geometrical 
structure  of  the  system  domain)  for  the  heat  equation  in  two  or  three  space 
dimensions.  These  so-called  domain  identification  problems  thus  constitute  an 
important  class  of  parabolic  system  identification  problems  related  to  space 
structures. 

In  this  paper  we  present  a  general  convergence/stability  (continuous 
dependence  on  observations)  framework  for  approximation  methods  to  treat 
parameter  identification  problems  (see  [BCK.])  involving  distributed  parameter 
systems.  This  new  parameter  estimation  convergence  framework  combines  a 
weak  version  of  the  system  (in  terms  of  sesquilinear  forms  in  the  spirit  used 
in  [BK1])  with  the  resolvent  convergence  form  of  the  Trotter-Kato 
approximation  theorem  [P],  [BK2J.  The  very  general  convergence  results 
depend  (in  addition,  of  course,  to  the  usual  properties  for  the  associated 
approximating  subspaces)  on  three  properties  of  the  parameter  (q  €  Q) 
dependent  sesquilinear  form  o( q)(-,  )  describing  the  system:  (A)  continuity 
(with  respect  to  the  parameter);  (B)  uniform  (in  the  parameter)  coercivitv:  and 
(C)  uniform  (in  the  parameter)  boundedness.  The  approach  permits  one  to 


give  convergence  and  stability  arguments  in  inverse  problems  under  extremely 
weak  compactness  assumptions  on  the  admissible  parameter  spaces  Q 
(equivalent  to  those  in  typical  variational  or  weak  approaches— see  (Bl), 
(BCR])  without  requiring  knowledge  of  smoothness  of  solutions  usually  a  part 
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of  the  variational  and  general  finite  element  type  arguments.  Thus  this 
approach  combines  in  a  single  framework  the  best  features  of  a  semigroup 
approximation  approach  (i.e.,  the  Trotter-fCato  theorem)  with  the  best  features 
of  a  variational  approach  (weak  assumptions  on  Q).  The  weakening  of  the 
compactness  criteria  on  Q  is  of  great  computational  importance  since  the 
constraints  associated  with  these  criteria  should  be  implemented  in  problems 
(see  [BI]  for  further  discussion  and  examples). 

An  additional  significant  feature  of  the  approach  described  in  this  paper 
is  that  for  first  order  systems  it  yields  directly  a  stronger  convergence  (in  the 
spatial  coordinate)  than  one  readily  obtains  without  extra  effort  in  either  the 
convergence  arguments  using  the  usual  operator  convergence  form  of  the 
Trotter-Kato  formulation  (e.g.,  see  Theorem  2.2  of  [BCK.]  and  Theorem  2.4  of 
[KLW])  or  the  finite  element  type  variational  arguments  (e.g.,  see.  Theorem  4.1 
of  [BR1]  and  the  related  remarks  in  [BK.L1]). 

While  the  approach  we  present  here  does  involve  coercivity  of  certain 
sesquilinear  forms  associated  with  the  system  dynamics,  its  applicability  is  nfit 
restricted  to  parabolic  systems  which  generate  analytic  semigroups.  As  we 
shall  demonstrate  in  Section  3  below,  it  can  be  used  to  treat  problems  in  which 
the  underlying  semigroup  is  not  analytic  (e.g.,  those  involving  Euler-Bernoulli 
equations  for  beams  with  Various  types  of  damping-viscous,  Kelvin-Voigt, 
spatial  hysteresis),  improving  substantially  on  some  of  the  currently  known 
results  for  these  problems.  With  appropriate  modifications,  this  theoretical 
framework  can  be  generalized  to  also  allow  an  elegant  treatment  of  problems 
involving  functional  partial  differential  equations,  e.g.,  beams  with  Boltzmann 
damping  (i.e.,  time  hysteresis)  --  see  [BFW], 
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2.  First  Order  Systems 


We  consider  first  order  systems  dependent  on  parameters  q  €  Q  described 
by  an  abstract  equation 


u(t)  =  A(q)u(t)  +  F(t,q) 
u(0)  =  u_(q) 


in  a  Hilbert  space  H.  The  admissible  parameter  space  Q  is  a  metric  space 
with  metric  d  and  for  q  €  Q,  we  assume  that  A(q)  is  the  infinitesimal 
generator  of  a  C0  semigroup  T(t;q)  on  H.  We  assume  that  we  are  given 
observations  u{  €  H  for  the  mild  solution  values  u(tjtq)  of  (2.1);  i.e.,  we  solve 
(2.1)  in  the  sense 


u(t;q)  =  T(t;q)u0(q)  +  T(t^;q)F(s,q)ds  (2.2) 

j0 

in  H.  We  then  consider  the  least  squares  identification  (ID)  problem  of 
minimizing  over  q  €  Q  the  functional 


J(q)  =  E  |u(tj;q)  -  Uj|2. 


Such  problems  are,  in  general,  infinite  dimensional  in  both  the  state  u  and 
the  parameter  q  and  thus  one  must  consider  a  sequence  of  computationally 
tractable  approximating  problems.  These  can  be,  for  our  purposes,  best 
described  in  terms  of  parameter  dependent  sesquilinear  forms  o(  q)(-,-) 
associated  with  (2.1)  or  (2.2)  (i.e.,  forms  which  define  the  operators  A(q)  in 
(2.1)).  For  more  details  and  examples  in  the  parabolic  case,  we  refer  the 
reader  to  (BK.1J.  Briefly,  let  V  and  H  be  Hilbert  spaces  with  V  continuously 


and  densely  imbedded  in  H.  Denote  a  family  of  parameter  dependent 
sesquilinear  forms  by  c<q):  V  x  V  I,  q  €  Q.  We  assume  that  o  possesses  the 
following  properties: 


(A)  Continuity:  For  q,q  €  Q,  we  have  for  all  4>,<l>  €  V 


l°(q)(4>,<M  “  o(q)(<M)l  «  d(q,q)|flv|fl. 


(B)  Coercivitv:  There  exist  cx  >  0  and  some  real  V0  such  that  for  q  €  Q, 


<p  €  V  we  have 


ofaXM)  +  VJ#*,  *  Cjlfl*. 


5.:  There  exist  c,  >  0  such  that  for  q  e  Q,  faip  e  V  we  have 


lo(q)(mi  (  c2mvmv. 


Under  these  assumptions,  a  defines  in  the  usual  manner  (e.g.,  see  [K],  [S]) 
operators  A(q)  such  that  o(q)(0,<M  ■  <-A(q)0,<J»H  for  0  €  dom(A(q)),  <p  €  V  with 
dom(A(q))  dense  in  V.  Furthermore,  A(q)  is  the  generator  of  an  analytic 
semigroup  T(t;q)  on  H  (indeed,  A(q)  is  sectorial  with  (VI  -  A(q))dom(A(q))  =  H) 
and  mild  solutions  of  (2.1)  possess  additional  regularity  (e.g.  see  Chap.  4  of  [P]  and 
Chap.  4  of  [5])  under  appropriate  assumptions  on  F.  Property  (B)  guarantees  that 
for  V  >  XQ  the  resolvent  operator  Rx(A(q))  =  (VI  -  A(q))'1  exists  as  a  bounded 
operator  on  H;  in  fact,  it  follows  readily  from  (B)  that  (assuming  ||H  S  k| -|v)  we 


Cl|Rx(A(q))<  «  <0, Rx(A(q))</»>H 
«  |<flH|Rx(A(q))«HH 
«  k*l<Mv|R  x( A(q))«Mv 
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and  hence  for  <|>  €  V 


while  for  4>  €  H  we  find 


|Rx(A(q))<Hv  <  —\Hy 

ci 


lRx(A(q))<MH  *  —  |#H. 

ci 


In  addition,  one  can  use  (B)  and  (A)  to  argue  that  q  -  R^(A(q))  is  a 
continuous  mapping  from  Q  to  V.  It  is  these  ideas  that  can  be  modified  to 
give  resolvent  convergence  in  the  approximation  schemes  which  we  introduce 


We  consider  Galerkin  type  approximations  in  the  context  of  sesquilinear 
forms  (e.g.,  see  [BK1]  for  further  details).  Let  HN  be  a  family  of  finite 
dimensional  subspaces  of  H  satisfying  PNz  -*  z  for  z  €  H  where  PN  is  the 
orthogonal  projection  of  H  onto  HN.  We  further  assume  that  HN  C  V  and  that 
the  family  possesses  certain  V-approximation  properties  to  be  specified  below. 
If  we  now  consider  the  restriction  of  o(q)(-,-)  to  HN  x  hn,  we  obtain 
operators  AN(q):  HN  •*  HN  which,  because  of  (B),  satisfy  a  uniform  dissipative 
inequality  and  can  be  shown  to  generate  semigroups  TN(t;q)  in  HN.  These  are 
then  used  to  define  approximating  systems  for  (2.2): 

uN(t;q)  =  TN(t;q)PNu0(q)  +  [  TN(t-s;q)PNF(s,q)ds.  (2.6) 

One  thus  obtains  a  sequence  of  approximating  ID  problems  consisting  of 


minimizing  over  Q 


JN(q)  =  E  |uN(t;;q)  -  us|2. 


1. 
ass 


vV  v.v?Xv :r.vV.v  " 


In  problems  where  Q  is  infinite  dimensional  (the  usual  case  in  many  inverse 
problems  of  interest),  on  must  also  make  approximations  QM  for  Q  (see  [BD] 
[BR2]  for  details).  To  include  this  aspect  of  the  problem  in  our  discussion 
here  does  not  entail  any  essential  mathematical  difficulties.  Since  it  would 
increase  the  notational  clutter  and  provide  no  new  mathematical  insight,  we 
do  not  pursue  the  theory  in  that  generality. 

To  obtain  convergence  and  continuous  dependence  (of  parameter  estimates 
with  respect  to  observations)  results  for  the  solutions  qN  of  minimizing  JN  in 
(2.7),  it  suffices  under  the  assumption  that  (Q,d)  is  a  compact  space  (see  [Bl]) 
to  argue:  for  arbitrary  (qN)  C  Q  with  qN  -•  q  we  have  uN(t;qN)  -*  u(t;q)  for  each 
t.  Under  reasonable  assumptions  on  F  and  uQ,.  this  can  be  argued  if  one  first 
shows  that  TN(t;qN)PNz  -*  T(t;q)z  for  arbitrary  qN  -*  q  and  z  e  H.  To  do  this 
one  can  use  several  versions  of  the  Trotter-  Kato  theorem  [P],  [BK2].  We  state 
precisely  the  "resolvent  convergence  form*  of  this  theorem  in  a  form  that  is 
general  enough  for  us  to  use  in  several  contexts  subsequently  in  this  paper. 

Let  X  and  XN,  N  =  1,2 .  be  Hilbert  spaces,  XN  C  X,  and  let  PN:  X  -  XN 

be  the  orthogonal  projection  of  X  onto  XN.  We  assume  the  XN,s  approximate 
X  in  the  sense  that  PNx  -»  x  for  all  x  €  X.  Then  we  have: 


Theorem  2.1:  Let  /1N  and  A  be  infinitesimal  generators  of  C0  semigroups  SN(t) 
and  S(t)  on  XN  and  X,  respectively ,  satisfying: 

(i)  There  exist  constants  w,  M  such  that  |SN(t)|  i  Me^*  for  each  N; 

(ii)  There  exists  X  €  p(A)  <\=1  p(/4N)  with  Re  X  >  w  such  that 
Rx(/tN)PNx  -*  Rx(>4)x  for  each  x  €  X. 
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(iii)  For  each  x  €  X,  SN(t)/>Nx  -» S(t)x  uniformly  in  t  on  any  compact  interval 
[O.tJ. 

Wc  may  assume  without  loss  of  generality  that  the  constants  u  and  M  of 
(i)  are  chosen  so  that  S(t)  also  satisfies  the  bound  in  (i). 

There  is  an  alternate  version  (we  shall  for  obvious  reasons  refer  to  this 
version  as  the  "operator  convergence  form")  of  the  Trotter-Kato  theorem  which 
has  been  frequently  used  [BCK],  [BC1},  [BC2],  [R]  in  parameter  estimation 
problems.  This  version  replaces  condition  (ii)  above  by  the  condition: 

(ii’)  There  exists  a  set  D  dense  in  X  such  that  for  some  X,  (XI  -  A) D  is  dense 
in  X  and  /tN/>Nx  -  Ax  for  all  xel. 

As  we  have  already  noted,  it  is  the  resolvent  convergence  form  of  this 
theorem  which  we  shall  wish  to  make  use  of  in  our  theoretical  framework 
since,  as  we  shall  show,  for  our  first  order  systems  condition  (ii)  will  follow 
readily  from  (A),  (B),  (C)  and  a  condition  on  how  HN  approximates  H.  The 
entire  theory  will  require  only  that  Q  be  compact  in  the  metric  d  used  in  the 
continuity  statement  (A)  for  the  sesquilinear  form.  This  will,  in  general,  be 
much  weaker  than  that  compactness  required  for  use  in  proofs  employing 
condition  (ii’),  since  the  requirements  on  "qN  -  q”  to  insure  AN(qN)PNz  ■*  A(q)z 
typically  involve  convergence  of  some  derivatives  of  the  qN.  For  example,  in 
first  order  parabolic  equations  containing  variable  coefficients  to  be 
estimated,  this  results  in  a  requirement  that  Q  be  compact  in  H1  whereas  use 
of  the  sesquilinear  form  approach  developed  in  this  paper  requires  only 
compactness  of  Q  in  C  or  L“  (see  [B1J  for  further  discussion  of  this  point). 
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Before  turning  to  our  convergence  arguments,  we  state  the  convergence 
properties  required  of  the  approximating  subspaces  HN  alluded  to  above. 
Throughout  our  discussions,  we  shall  assume: 

(Cl)  For  each  z  e  V,  there  exists  zN  €  HN  such  that  |z  -  zN|v  -»  0  as 
N  -  “ 

Theorem  2.2:  Let  conditions  (A),  (B),  (C)  and  (Cl)  hold  and  qN  -*  q  in  Q.  Then 
for  X  =  XQ,  Ry(AN(qN))PNz  -*  Ry(A(q))z  in  the  V  norm  for  any  z  €  H. 

Proof:  First  note  that  since  AN(qN)  results  from  the  restriction  of  o  to  HN  x 
Hn,  we  can  choose  X  -  XQ  in  (B)  so  that  X  e  p(A(q))  n£=1  p(AN(qN))  with  X  >  w 
(and  indeed  bounds  similar  to  those  of  (2.4),  (2.5)  hold  for  Ry(AN(qN))PN). 

Let  z  €  H  be  arbitrary  and  for  notational  convenience  put  w  =  w(q)  = 
R^(A(q))z  and  wN  =  wN(qN)  s  Ry(AN(qN))PNz.  Since  w  e  dom(A(q))  C  V,  we 
may  use  condition  (Cl)  to  define  a  sequence  (wN)  satisfying  |wN  -  w|v  -  0  as 


We  wish  to  show  |wN  -  w|y  -  0.  If  we  argue  that  |wN  —  wN|y  •*  0,  the 
desired  results  follow  immediately  from  the  triangle  inequality. 


Let  zN  =  wN  -  wN  so  that  zN  €  Hw  C  V.  Using  <•,->  to  denote  the  inner 
product  in  H,  we  find 


o(q)(w,zN)  =  <— A(q)w,zN>  =  <(X  —  A(q))w,zN>  —  X<w,zN> 

=  <z,zN>  -  X<w,zN> 


o(qN)(wN,zN)  -  <— AN(qN)wN,zN>  =  <(X  -  AN(qN))wN,zN>  -  X<wN,zN> 

=  <PNZ,ZN>  -  X<WN,ZN> 

=  <z,zN>  -  X<wN,zN>. 

Hcncc  we  have 

o(q)(w,zN)  =  o<qN)(wN,zN)  +  X<wN  -w,zN>. 

Using  this  with  (B)  we  obtain 

cjz1^  S  o(qN)(zN,zK)  +  X|zN|2  =  o(qN)(wN,zN)  -  o(qN)(wN,zN)  +  X|zN|j 
=  o(q)(w,zN)  +  X<w-wN,zN>  ~  o( qN)(wN,zN)  +  X|zN|^ 

=  o(q)(w,zN)  -  c<qN)(w,zN)  +  o(qN)(w-wN,zN) 

+  X(<w-wN,zN>  +  |zN|^). 

If  we  use  (A)  and  (C),  we  thus  have 

ct|zN|2  S  d(q,qN)|w|v|zN|v  +  c2|w-wN|v|zN|v  +  X(<w-wN,zN>  +  |zN|2). 

Finally,  since 

<w-wN,zN>  =  <w— wN,zN>  +  <wN— wN,zN>  =  <w-wN,zN>  —  <zN,zN>, 

the  above  inequality  may  be  written 

c,|zNl^  S  d(q,qN)|w|v|zN|v  +  c2|w  -  wN|v|zN|v  +  X<w  -  wN,zN> 

«  d(q,qN)|w|v|zN|v  +  c2|w  -  wN|v(zN|v  +  |X||w  -  wN|H|zN|H. 

Thus,  using  |-|H  S  k|  |v  we  find 


c^ly  «  d(q,qN)|w|y  +  (c2  +  |X|k2)|w  -  wN|v, 


which,  from  the  definition  of  {wN}  and  the  fact  that  qN  -*  q,  yields  the 
desired  convergence  zN  =  wN  —  wN  -•  0  in  V  norm. 

Since  V  norm  convergence  is  stronger  than  H-norm  convergence,  we  can 
use  the  results  of  this  theorem  along  with  Theorem  2.1  for  X  =  H,  XN  =  HN, 
pN  _  pN^  _  AN(qN),  A  _  A(q)  to  obtain  immediately  that  TN(t;qN)PNz  - 
T(t;q)z  for  z  G  H,  the  convergence  being  in  the  H  norm.  This  along  with  (2.2) 
and  (2.6)  can  be  used  to  argue  the  desired  convergence  uN(t;qN)  -  u(t;q)  in  the 
H  norm  (assuming,  of  course,  appropriate  smoothness  of  F  and  u„  in  q  e  Q). 
However,  with  no  additional  assumptions  and  only  a  little  extra  effort,  one 
can  obtain  this  convergence  in  the  V  norm.  As  we  shall  explain  below,  this 
stronger  convergence  is  often  immensely  useful  in  parameter  estimation 
problems  where  the  observations  may  be  continuous  in  the  V  norm,  but  not  in 
the  H  norm. 

Theorem  2.3:  Under  the  hypotheses  of  Theorem  2.2,  we  have  for  each  z  €  H, 
TN(t;qN)PNz  -  T(t;q)z  in  the  V  norm  for  t  >  0,  uniformly  in  t  on  compact 
subintervals. 

Proof:  To  establish  these  results,  we  first  apply  the  Trotter-Kato  theorem  in 
the  space  X  =  V  with  XN  =  HN.  Let  P^:  V  -»  HN  denote  the  orthogonal 
projection  of  V  onto  HN  in  the  V  inner  product.  Then  hypothesis  (Cl)  implies 
Pyf  -  f  in  the  V  norm  for  any  f  G  V  as  well  as  PNf  -  f  in  the  H  norm  for 
any  f  G  H. 

To  carry  out  our  arguments,  we  shall  employ  several  bounds  for  resolvents 
similar  to  but  sharper  than  those  of  (2.4),  (2.5).  These  bounds  follow  from 


,«<  •  •*»  '„<«*! 


results  in  Tanabe  (c.f.  Lemma  6.1  of  Chapter  3  in  [T]);  we  note  that  our 
operator  A  can  be  extended  to  A:  V  -*  V*  as  in  [T]  so  that  the  results  given 
there  are  applicable  as  used  below.  (In  [T],  the  condition  (B)  is  assumed  to 
hold  with  XQ  =  0  and  we,  without  loss  of  generality,  state  our  bounds  for  this 
case.)  The  first  bound  we  shall  use  yields  for  f  €  V 

|Rv(AN(qN»f|v  <  j^jT7i-in„  (2  8) 

where  c  is  independent  of  N.  Using  this  with  f  =  P^z  -  PNz  for  z  e  V  along 
with  the  results  of  Theorem  2.2  we  obtain  Rx(AN(qN))PyZ  -*  Rx(A(q))z  in  V 
norm  for  any  z  €  V.  This  is  (ii)  of  Theorem  2.1  with  X  =  V,  XN  =  HN  and 
PN  =  Py. 

To  obtain  the  uniform  stability  bound  (i)  of  Theorem  2.1  in  the  V-norm, 
we  make  use  of  another  bound  that  follows  from  [T],  From  the  last  bound  of 
Lemma  6.1,  Chap.  3  of  [T]  one  may  readily  argue  that  for  f  €  V 


as? 


|Rl(AN(qN))f|v  <  -HC|V 


where  Mx  is  independent  of  N.  The  arguments  behind  (2.9)  involve  using 
(3.51)  of  [T]  for  A*  and  observing  that  for  7  €  V,  v  e  V*  arbitrary  so  that 
RX(A)7  €  V  =  V**  we  have  (using  the  usual  notation  for  the  duality  product) 


(Rx(A)7,v)v..  v.  =  (v,Rx(A)7)v,  y  =  (R X(A*)v,7)v.  v 

M. 

«  |Rx(A*)v|v.i7iv  s  -H7IvMv- 


v.v> 
*►  -  -  .% 


& 


Then  |Rx(A)7lv.,  «  |7lv  and  hencei  identifying  V  and  V**,  we  obtain 
M, 

|RX(A)7IV  <  y^-lTlv  for  7  €  v- 


- 1 J- 


We  may  use  the  bound  (2.9)  in  the  identity 


TN(t;qN)  =  -V  [  euRx(AN(qN))dX, 
2m  Jr 


(2.10) 


where  r  is  a  contour  (about  X.)  similar  to  that  given  in  [P,  p.  63],  to  argue  for 


€HN 


|TN(t;qN)z|v  S  McV°‘|zlv 


(2.11) 


where  M  is  independent  of  N.  Thus  one  obtains  |TN(t;qN)|v  S  Me^o*.  (These 
results  are  essentially  given  in  [T];  one  just  needs  to  check  the  arguments  to 
ascertain  that  the  constant  M  depends  on  AN  only  through  the  bounds  ct  and  c2 
of  conditions  (B)  and  (C)  and  the  sector  angle  6  for  AN  --  where  6  again 
depends  only  on  c2  and  cr) 

We  may  then  apply  Theorem  2.1  as  indicated  to  obtain  TN(t;qN)Pyz  •* 
T(t;q)z  in  V  norm  for  every  z  €  V.  It  remains  to  argue  that  TN(t;qN)PNz  - 
T(t;q)z  in  V  norm  for  every  z  €  H. 

If  we  use  the  bound  (2.8)  in  (2.10),  arguments  similar  to  those  used  to 
establish  (2.11)  can  be  made  to  give 


|TN(t;qN)z|v  «  Me  °V1/2|z|H,  t  >  0,  z  €  HN  , 


(2.12) 


where  again  M  is  independent  of  N. 

Recalling  that  V  is  dense  in  H,  then  given  z  G  H  and  t  >  0,  we  may 
choose  zv  €  V  such  that  |zy  -  z|H  <  e.  Then  noting  that  T(t;q)z  e  V  for  t  >  0, 
wc  have 
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|TN(t;qN)PNz  -  T(t;q)z|v  <  |T"(t;q")(P"z  -  P$zv)|v 

+  |TN(t;qN)PyZv  -  T(t;q)zv|v  +  |T(t;q)(zv  -  z)|v. 

From  (2.12)  we  see  that  the  first  term  in  the  right  side  of  this  inequality  is 
bounded  by  Mcxott'I/2|PNz  -  Pvz\4h  which,  as  N  -  ®  approaches 

Mexott'1/,2|z-zvlH.  The  third  term  is  also  bounded  by  a  similar  expression, 
while  the  second  term  approaches  zero  as  N  -»  08  by  our  previous  arguments. 
We  thus  have  established  the  convergence  claimed  in  the  statement  of  Theorem 


N/».„NwpN_  _  PN, 


Among  the  examples  that  can  be  treated  immediately  with  the  above 
theory  are  the  usual  parabolic  systems  (see  [BK.1],  [L]).  To  illustrate  these 
ideas,  consider  the  estimation  problem  for  the  standard  Dirichlet  boundary 
value  problem  for  one-dimensional  parabolic  systems.  That  is, 

9u  3  f  6u  1  3 

iT~  ST  Lq*  srj  *  +  ^  on  °-(<U) 

u(t,0)  =  u(t,  1 )  =  0, 


with  q  =  (q1,q2,q3)  to  be  chosen  (via  a  least  squares  criterion)  from  Q,  a 
compact  subset  of  C(n)  x  C(H)  x  C(fl).  The  state  spaces  are  H  =  H0(D),  V  = 
Hj(fi)  and  the  weak  form  of  the  equation  is  given  by 


<z(t),0>  +  o(q)(z(t),</>)  =  0,  ijie  V 


where  the  sesquilinear  form  is  defined  as 


=  <qjD9,Dii»>  +  <q2<f,D0>  -  <q3<p,0>. 


with  D  =  9/3x.  One  can  readily  verify  that  conditions  (A)  and  (C)  hold, 
while  the  coercivity  condition  (B)  is  valid  if  we  assume  elements  of  Q  satisfy 
q2(x)  >  v  >  0  for  some  constant  v. 

In  this  case  one  can  specify  the  domain  of  A(q)  by 

dom  A(q)  =  [<p  €  H*(Q)  |  qjDtp  +  q2<p  €  H1^)} 

and  A(q)(p  =  D(q1D<p)  +  D(q2<f>)  +  q3<p.  For  approximating  elements  one  can  use 
either  piecewise  cubic  or  piecewise  linear  B-splines  modified  to  satisfy  the 
Dirichlet  boundary  conditions. 

We  note  that  the  theory  above  only  requires  compactness  of  Q  in  the 
[C(fi)]3  topology  along  with  relatively  weak  smoothness  assumptions  when 
compared  with  other  approaches  (see  [BK3],  [BKL1],  [BKL2],  [BCK]). 
Furthermore  we  obtain  V  =  Hj(Q)  convergence  of  the  states.  Hence  the  theory 
is  complete  for  least  squares  criterion  involving  pointwise  (in  x)  observations 
of  the  state.  This  is  obtained  for  little  extra  effort  when  compared  with  the 
efforts  usually  required  to  obtain  the  stronger  convergence  (see  [BCK.], 
[BKL1]). 

Other  parabolic  examples  of  great  practical  importance  can  be  readily 
treated  in  the  context  of  this  framework.  For  example,  the  2-D  thermal 
tomography  problem  mentioned  in  the  Introduction  and  discussed  in  [BKO] 
results  in  a  boundary  identification  problem  for  parabolic  systems.  With  a 
standard  transformation  these  problems  can  be  readily  treated  theoretically  and 
computationally  with  the  approach  of  this  paper.  Details  can  be  found  in 
[BKO]. 

Another  class  of  problems  of  interest  involves  estimation  of  coefficients 


in  the  Fokker-Planck  or  forward  Kolmogorov  equations  (see  [B2])  for  the  case 
where  the  Markov  transition  process  for  growth  is  time  invariant  in  size/age 
structured  population  models.  In  this  case  the  equations  have  the  form 

3u  3  32 

iT  +  3^QlU)  =  dJ (q’u)  "  q3u’  x°  <  x  <  xv 
with  boundary  conditions 


f  3  "I x— xo  fxi 

(q2u)J  =  j  q4(Ou(t,Ud( 

xo 

r  3  ix=xi 

LqiU  ~3T(q2U)J  =0‘ 


The  associated  sesquilinear  form  is  given  by  (here  H  =  H^x^Xj),  V  = 

H^Xq.Xj)) 

o(q)(<P,<W  =  “<q1<P  “  D(q2<p),Dil»>  -  <Kx0)R(q4)(<p)  +  <qs‘P»^> 

with  R(q4)(<p)  =  f*1  q4(5)<p(OcU.  Under  appropriate  assumptions  on  Q  C  L^n) 

*o 

*  *  L<a(n)  x  Ljn),  n  =  (x^x^.  (see  [B2]  for  details)  one  can  readily  argue 

that  hypotheses  (A),  (B),  (C)  hold.  Thus  these  problems  also  fall  within  the 
purview  of  the  theory  developed  here. 
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3.  Second  Order  Systems 


The  ideas  developed  in  the  previous  section  can  be  applied  to  parameter 
estimation  problems  involving  second  order  systems  of  the  form 


u(t)  +  B(q)u(t)  +  A(q)u(t)  =  f(t) 


in  a  Hilbert  space  H  where  the  operators  A(q)  and  B(q)  are  defined  via 
parameter  dependent  sesquilinear  forms  in  a  manner  similar  to  that  of  Section 
2.  We  use  the  general  approach  given  in  [S];  we  again  assume  we  are  given  a 
Hilbert  space  V  C  H  that  is  continuously  and  densely  imbedded  in  H.  Let 
ox(q):  V  x  V  -»  <E  be  a  symmetric  sesquilinear  form  satisfying  conditions  (A) 
and  (C)  of  Section  2  and  the  V  coercivity  condition  (B)  with  XQ  =  0.  Then  for 
each  q  €  Q  we  may  define  continuous  linear  operators  A(q):  V  -*  V*  given  by 


®j(qXM)  -  (A(q)<M)v.  v 


for  <P,4>  €  V,  where  as  customary,  we  have  identified  the  pivot  space  H  with  its 
dual  H*.  i.e.,  V  C  H  =  H*  C  V*  and  the  duality  product  (  -,  -)v.  v  is  the  unique 
extension  by  continuity  of  the  scalar  product  <■,•>  of  H  from  H  x  V  to 
V*  x  V  (e.g.,  see  [S]  for  details).  Of  course,  as  in  Section  2,  we  may  also  view 
A(q)  as  a  densely  defined  operator  in  H  where  at( q)(4>,4>)  =  <A(q )<P,<P>H  for 
<p  €  dom  A(q),  <p  €  V.  The  sesquilinear  form  Oj  and  its  associated  operator  A 
will  correspond  to  the  "stiffness"  operator  in  the  examples  we  treat  below. 

A  second  sesquilinear  form  will  give  rise  to  the  "damping"  term  in  our 
examples.  We  assume  we  arc  given  a  (not  necessarily  symmetric)  sesquilinear 
form  o,(q):  V  x  V  -»  C  satisfying  conditions  (A)  and  (C)  and  the  semicoercivitv 


/•.■Vi 

v.  .•> 


condition 


(B1)  H-semicoercivitv:  There  exists  b  >  0  such  that  for  q€Q  and  <l>  €  V 
we  have 

o2(q)(W)  >  (3.3) 

As  before,  under  condition  (C)  on  aJt  we  may  define  a  continuous  linear 
operator  B(q):  V  -»  V*  given  by 

°j(q)($,tM  =  (B(q)0»v.  v.  (3.4) 

Again,  we  may  alternatively  view  B(q)  as  being  densely  defined  in  H. 
However,  since  we  may  wish  in  general  to  view  equation  (3.1)  as  an  equation 
in  V*  (again  following  classical  formulations  [S],  [L],  [T])  we  shall  interpret 
A(q)  and  B(q)  as  members  of  2(V,V*).  We  note  that  from  conditions  (B)  and 
(C)  on  Oj,  the  form  Oj(q)  is,  for  each  q  €  Q,  equivalent  to  the  norm  inner 
product  in  V;  indeed,  we  can  use  o2(q)  to  define  a  parameter  dependent 
equivalent  inner  product  in  V.  We  shall  use  to  denote  the  Hilbert  space 
consisting  of  the  elements  of  V  equipped  with  this  inner  product. 

We  shall  rewrite  equation  (3.1)  in  first  order  form  and  to  that  end  we 
define  the  product  spaces  V  =  V  x  V  and  K  =  V  x  H.  If  <-,->v  denotes  the 
inner  product  in  V,  we  may  use  and  o2  to  define  a  sesquilincar  form 
a  :  V  x  V  **  <E  given  by 

o(q)((u,v),($,<(0)  =  -<v,$>v  +  Oj(q)(u,0)  +  o2(q)(v,^).  (3.5) 

We  may  then  write  equation  (3.1)  in  weak  or  variational  form  as 

<w(t),X>K  +  <*q)(w('),X)  =  <F(t),X>K  •  X  G  V,  (3.6) 


where  w(t)  =  (u(t),v(t)),  x  =  (0,0),  F(t)  =  (f(t),0)  and  <-,•  >#  is  the  usual  product 
space  inner  product.  Alternatively,  we  may  write  equation  (3.1)  as 


w(t)  =  ^(q)w(t)  +  F(t) 


where  A( q)  is  the  operator  associated  with  the  sesquilinear  form  a  in  the  usual 
manner. 

Following  this  approach,  we  define  in  If  =  V  x  H  the  operator 


^(q)  = 


-A(q)  -B(q) 


on  dom  >4(q)  =  {(0,0)  €  K  |  0  €  V  and  A(q)0  +  B(q)0  €  H)  C  V.  We  can  readily  argue 
that  dom  /4(q)  is  dense  in  If  =  V  x  H  and  that,  under  the  conditions  on  ot  and  o2,  for 
X  >  0  the  range  of  XI  -  A( q)  is  X  (see  the  arguments  given  below).  Moreover,  A(q) 
is  dissipative  in  Xq  h  Vq  x  H  since  for  (0,0)  €  dom  A(q) 

</t(q)(0,0),(0,0)>x  =  o1(q)(0,0)  -  Oj(q)(0,0)  -  o2( q)(0,0) 

<  -  b(01*  (  0. 

Thus  by  the  Lumer-Phillips  theorem  (e.g.  see  Chap.  1,  Thms.  4.3,  4.6  of  [P])  we  find 
that  /4(q)  generates  a  semigroup  of  contractions  on  Vq  x  H.  Since  Vq  and  V  possess 
equivalent  (uniformly  in  q  e  Q)  norms,  we  have  that  /4(q)  is  the  infinitesimal 
generator  of  a  CQ-semigroup  T(t;q)  on  X  =  V  x  H.  If  the  form  o2  satisfies  (3.3)  with 
b  >  0,  then  one  can  show  that  this  semigroup  is  uniformly  exponentially  stable, 
i.c.,  |T(t;q)j  <  Mc'u)t  for  some  M  ?  1  and  w  >  0.  Furthermore,  if  the  H 
scmicoercivity  of  (3.3)  is  replaced  by  V  coercivity,  i.e.,  b  >  0  and  the  H  norm  is 
replaced  by  the  V  norm  in  (3.3),  then  T(t;q)  is  an  analytic  semigroup  on  X  =  V  «H. 


s 

-  J'** 

•  «<  *  ■ 

_  •  _  »  _ 


;s 


vVv 

VV 

•>> 


>' v  . 

JP>« 


Since  we  have  o(q)(x,0  =  ~C4(q)x,(  >jf  for  x  €  dom  A( q)  and  5  e  V  = 
V  x  V,  one  might  be  tempted  to  apply  the  arguments  and  results  of  Section  2 
to  this  sesquilinear  form  with  V  and  K  playing  the  role  of  the  spaces  V  and  H 
in  the  theory  of  Section  2  (i.e.,  treat  this  system  generated  by  <4(q)  as  a  first 
order  system  in  K  and  apply  without  modification  the  theory  in  Section  2).  If 
one  did  this,  then  conditions  (A),  (B)  and  (C)  must  be  satisfied  by  o(q)  using 
the  norms  of  K  and  V.  Conditions  (A)  and  (C)  pose  no  difficulty  under  the 
assumptions  of  (A)  and  (C)  on  ox  and  o2  given  above.  However  to  argue  that 
o(q)  is  V  coercive  (condition  (B)  of  Section  2),  one  finds  that  the  H 
semicoercivity  condition  (3.3)  on  o2(q)  must  be  strengthened  to  V  coercivity. 
While  some  damping  forms  of  interest  (e.g.  strong  Kelvin-Voigt  damping  --  we 
discuss  this  below)  do  satisfy  this  stronger  condition,  several  important  types 
of  damping  lead  to  sesquilinear  forms  that  don’t.  Hence  we  shall  modify  the 
arguments  behind  Theorem  2.2  in  order  to  enable  us  to  treat  the  more  general 


We  shall  throughout  our  discussions  henceforth  assume  that  at  satisfies 
conditions  (A),  (B)  with  XQ  =  0,  and  (C),  while  o2  satisfies  conditions  (A),  (B1),  (C). 
Again  we  shall  be  interested  in  a  resolvent  convergence  form  of  the  Trotter-Kato 
theorem  and  hence  must  consider  the  resolvent  of  the  operator  /4(q)  given  in 


To  motivate  our  discussions,  consider  for  X  >  0  the  equation  in  K  given  by 
X  =  Rjk(/t(q))(  for  X  =  (9,0),  (  =  (n,7).  This  is  equivalent  to  solving  for  X  € 
dom(/4(q))  in  the  equation  (XI  -  A(q))x  =  (  which  may  be  written 


\<p  —  Ip  =  T) 

(3.9) 

\<l>  +  A(q)<p  +  B(q)<l>  =  y  . 

If  we  substitute  the  first  equation  <J>  =  X<p  -  n  of  this  system  into  the  second 
we  have 


x2<p  +  A(q)<p  +  XB(q)<?  =  y  +  xn  +  B(q)n 


(3.10) 


which  must  be  solved  for  <p  €  V.  This  suggests  that  we  define,  for  X  >  0,  the 
associated  sesquilinear  form  o^(q):  V  x  V  -*  C  given  by 

°x(q)(<p,0  =  ^<<P,5>h  +  0i(q)(<P.C)  +  Xo2(q)(<p,0.  (3.11) 

Since  al  is  V-coercive  and  o2  is  H-semicoercive,  we  find  for  <p  €  V 

0\(q)(<P,9)  =  +  OjfaX^.ip)  +  Xo2(q)(<p,(p) 

>  x2Mh  +  ciMy  +  xbltpl*. 


Hence  for  <p  €  V 


0x(q)(<i>.<p)  >  +  xm2  >  CjM2 


(3.12) 


where  X  =  (X2  +  Xb)  >  0.  Thus,  for  X  >  0,  the  form  is  V-coercivc  and  the 
equation  (3.10)  is  solvable  for  <p  e  V  for  any  given  £  =  (n,y)  in  X  =  V  x  H.  It  follows 
then  that  defining  <l>  using  the  first  equation  in  (3.9),  we  may,  for  any  {  =  (n,y), 
find  an  element  x  =  (<p,0)  in  dom(^(q))  which  solves  (3.9);  i.e.  for  X  >  0,  Rx(/!(q)) 
exists  as  an  element  in  2(H). 

The  coercive  inequality  (3.12)  will  be  the  basis  of  our  convergence 
arguments. 
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Wc  consider  next  a  Galerkin  type  approximation  scheme  for  equation  (3.7). 
As  in  Section  2,  we  assume  approximation  subspaces  HN  C  V  satisfying 
condition  (Cl).  Then  for  any  w  e  V  =  V  x  V  and  each  N,  there  exists  wN  e 
Hn  x  Hn  satisfying  as  N  -  • 

|wN  -  w| y  -  0.  (3.13) 

To  define  the  approximating  systems,  we  consider  o(q)  defined  in  (3.5) 
restricted  to  JfN  =  HN  x  HN  and  obtain  in  the  usual  manner  the  operators 
/lN(q):  Kn  -♦  Kn.  Since  (3.12)  holds,  we  have  immediately  that  Ry(^N(q))  exists 
in  X(Kn)  for  X  >  0.  Denoting  by  PN  the  orthogonal  projection  of  If  onto  JfN,  we 
may  then  prove  the  following  convergence  results  which  is  analogous  to  that 
of  Theorem  2.2  for  first  order  systems. 

Theorem  3.1.  Suppose  that  conditions  (A),  (B)  with  XQ  =  0,  (C)  for  or  (A),  (B'), 
(C)  for  o2  and  (Cl)  hold  and  let  qN  -*  q  in  Q.  Then  for  X  >  0  we  have 
RxUN(qN))PNt  **  R^(q))(  in  the  V  norm  for  any  5  €  X  . 

Proof.  Let  ?  €  Jf  be  arbitrary  and  put  w  =  w(q)  =  Rx(/4(q))(,  wN  =  wN(qN)  = 
R^(/4N(qN))PNL  Let  w  =  (<p»,  (  =  (n,7)  so  that  equation  (3.9)  is  satisfied  by 
this  pair.  Also,  letting  wN  =  PN(  =  (^N,7N),  we  have  that  this  pair 

satisfies  (3.9)  with  q  =  qN.  In  particular,  using  (3.10)  and  (3.11)  we  find  for 


5  €  V 


0x(q)(<P,O  =  <7,5  >H  +  x<7J,5>h  +  o2(q)(n,0 
ox(  qN)(<PN,0  =  <7N,5  >h  +  x<^N,5>h  +  ^(qN)(nN,5). 
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Since  w  =  w(q)  e  dom(/4(q))  c  V,  we  may  choose  wN  =  (^,0^)  in  kn 
satisfying  (3.13).  Thus,  choosing  5  =  £N  =  (pN  -  <pN  in  HN  C  V,  we  may  use  the 
above  equalities  to  obtain 


ox(q)(<P,SN)  "  °x(qN)(<pN,SN)  =  <7^N,;N>h  +  X<n-nN,QN>H 


+  o2(q)KSN)  -  o2(qN)(nN,;N). 


(3.14) 


We  use  this  to  argue  that  QN  -  0  in  V. 

Applying  the  estimate  (3.12)  to  the  element  =  <fN  -  <pN  and  then  using 
(3.14)  we  find 


cjt*^  +  X|SN&  S  ox(qN)(;N,5N) 

=  ox(qN)(9N,5N)  -  ox(q)(<p,5N)  +  ox(q)(<p,CN)  "  ox(qN)(9N,£N) 

=  <7N-7»5n>h  +  x<nN-n,5N>H  +  o2(q)(n,CN)  +  o2(qN)(nN,tN) 

+  ox(q)(9,5N)  -  ox( qN)(9N,5N)- 

The  last  four  terms  on  the  right  side  of  this  equality  may  be  written  as 


(3.15) 


o2( q)(nN-n,SN)  -  o2(q)(nN,5N)  +  o2(qN)(nN,;N) 

+  ox(q)(<p-<PN,;N)  +  °x(q)(?N,^N)  -  ox(qN)(;N,^N). 


(3.16) 


Using  the  definition  (3.11)  of  ox  and  the  boundedness  and  continuity  in  q  of 
Oj,  o2  (condition  (C)  and  (A))  along  with  |  |H  S  k|  |v  wc  may  bound  from  above 
the  expression  in  (3.16)  by 


c2|nN-n|vKN|v  +  d(q,qIN)|n^|viriv  +  C|«rVlvKwlv 


Nm„N,  ir  N, 


,r  N| 


+  (l  +  V)d(q,qN)|<PNlvi;N|v. 


(3.17) 


where  C  =  X2k2  +  c2  +  Xc2.  We  observe  that  since  PN  is  the  orthogonal 
projection  of  K  onto  KN,  hypothesis  (Cl)  implies  that  (oN,yN)  =  PN(n,y)  -*  (n,y) 
in  It  =  V  x  H.  Thus,  |nN|v  «  K.  for  some  constant  K.  Since  <pN  -•  (p  in  V,  we 

a  ij 

may  also  assume  |<p  |v  (  K. 

A 

Recalling  that  X  >  0,  we  can  combine  (3.15),  (3.16),  (3.17)  and  the  above 
observations  to  obtain  the  estimate 

CjKn|v  «  k|/SlH  +  (Xk  +  c2)|nN-n|v  +  K(2+X)d(q,qN)  +  C|<rvNlv, 

and  thus  £N  -*  0  in  V.  That  is,  <pN  -  <^N  -*  0  in  V  and  hence  <pN  -*  <p  in  V. 

Recalling  that  (<pN,ij^),  (nN,7N)  and  (%ip),  (n,y)  satisfy  (3.9)  with  qN  and  q, 
respectively,  we  have  that 

ip^  =  x/  -  nN  and  ip  =  \<p  -  n. 

Thus  we  find  immediately  that  ^  -  <p  in  V.  This  completes  the  proof  that 
wN  =  (vN,<pN)  ■*  w  »  (<p,ip)  in  the  V  norm. 

Let  TN(t;qN)  denote  the  CQ  semigroup  generated  by  /4N(qN)  in  ttN.  Then 
applying  Theorem  2.1  with  X  =  K  ,  XN  =  KN  and  /1N  =  /lN(qN),  we  obtain 
immediately  from  the  above  discussions  the  following  desired  convergence 
results. 

Theorem  3.2.  Under  the  hypotheses  of  Theorem  3.1,  we  have  for  each  (  €H  and 
t  >  0,  TN(t;qN)PN(  -  T(t;q)(  in  K,  with  the  convergence  being  uniform  in  t  on 
compact  subintervals. 

The  needed  convergence  results  for  parameter  estimation  problems  follows 
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directly  from  this  theorem.  In  particular,  for  solutions  of  (3.7)  and  their 
approximations  we  have  (assuming,  of  course,  appropriate  smoothness  in  q  of 
initial  data  and  the  perturbation  function  F) 


uN(t;qN)  -  u(t;q)  in  V  norm, 
vN(t;qN)  -  v(t;q)  in  H  norm. 


» 
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To  illustrate  the  types  of  problems  that  can  be  investigated  in  the  context  of 
the  theoretical  framework  developed  in  this  section,  we  consider  a  cantilevered 
Eulcr-Bcrnoulli  beam  of  unit  length  and  unit  linear  mass  density.  The  transverse 
vibrations  can  be  described  by  an  equation  of  the  form  (3.1)  with  the  operators 
A(q)  and  B(q)  chosen  appropriately  (see  e.g.  [BCR],  [BR2],  [R]  for  detailed 
equations).  We  assume  that  the  end  at  x  =  0  is  fixed,  the  end  atx  =  lisfrec.  Then 
the  state  spaces  can  be  chosen  as  H  =  H°(0,1)  and  V  =  H2(0,1)  =  (<p  e  H2(0,i)  |  <p(0)  = 
D<p(0)  =  0).  The  stiffness  sesquilinear  form  Oj  is  given  by 


o1(q)(<p»  =  f  EID 

where  we  assume  the  stiffness  coefficient  qt(x)  =  EI(x)  satisfies  q^x)  J  a  >  0 
and  is  in  Lm(0,l).  The  operator  A(q)  in  (3.1)  is  given  by  A(q)<p  =  D2(ElD2<+>) 
with 

dom  A(q)  =  (<p  €  V  |  EIDz<p  e  H2(0,1),  EIDV.l)  =  D(ElD2ip)(  1 )  =  0). 


If  we  denote  by  q2  the  damping  parameter  to  be  discussed  in  the  following 
examples,  then  the  admissible  parameter  set  Q  can  be  taken  as  a  compact 


subset  of 


V** 


Q  =  {q  =  (q^dj)  I  Q,  €  LJO,l),  qjx)  *  a  >  0). 


With  this  formulation,  a  number  of  important  damping  mechanisms  can  be 
readily  treated. 


Viscous  damping:  In  this  case  air  or  fluid  damping  is  usually  assumed 
proportional  to  the  velocity  of  displacement  so  that  the  term  B(q)u(t)  in  (3.1) 
has  the  form  bu(t)  for  q,(x)  =  b(x).  The  damping  sesquilinear  form  a.  is  given 


0j(q)(<P,<M  *  <btp,0  >H 


which,  for  b(x)  >  0,  satisfies  the  H-semicoercivity  condition  (B1)  but  of  course 
is  not  V -coercive.  The  domain  of  the  operator  in  (3.8)  is  given  by  dom  A( q)  = 


dom  A(q)  x  V. 


)ing:  For  these  models  the  damping  moment  is  postulated  as 


being  proportional  to  the  strain  rate  and  hence  the  damping  term  has  the  form 
D2(cDID2u(t))  where  q2(x)  =  cpI(x).  The  associated  sesquilinear  form  is  given 


o2(  q)(<P,<W  =  <cDIDVD2*>H 

which  for  cDI(x)  t  0  >  0  satisfies  a  V-coercive  condition  and  hence  of  course 
the  H-scmicoercivity  of  hypothesis  (B’).  Letting  M  =  EID2ip  +  cdID20  denote 
the  "moment"  of  the  beam,  we  have  the  domain  of  A(q)  given  by 

dom  -4(q)  =  {(<p,0)  €  V  x  H  |  <i>  €  V,  Me  H2(0,I),  M(l)  =  DM(I)  =  0). 


If  we  only  have  cDI(x)  ?  0,  then  only  condition  (B1)  is  satisfied  but  our 


$ 


K 
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theoretical  framework  is  still  applicable.  We  note  that  results  of  the 
framework  then  provide  a  distinct  improvement  over  the  results  for  the 
problems  given  in  [BCR]  (where  smoothness  assumptions  on  solutions  must  be 
hypothesized  as  well  as  q2  >  6)  and  [KG]  (where  a  stronger  and  less  readily 
characterized  topology  than  that  of  C(0,i)  or  L^O,!)  must  be  used  for  the 
compactness  assumption  on  Q). 

Spatial  hysteresis  damning:  This  damping  model,  which  has  been  suggested 
recently  and  investigated  by  Russell  in  [RU],  is  based  on  physical  foundations 
that  appear  particularly  appropriate  for  fiber  reinforced  composite  beams, 
e.g.,  beams  of  composite  materials  in  which  strands  of  fibers  are  lengthwise 
bound  in  an  epoxy  matrix.  The  damping  term  has  the  form  B(q)u(t)  = 
D((G(q)  -  vI)Du(t))  where  G(q)  is  a  compact  operator  in  H°(0,1)  defined  by  a 
symmetric  kernel  q2(x,y)  =  b(x,y)  =  b(y,x)  >0  a? 

(G<p)(x)  =  [  b(x,y)«p(y)dy,  <p  e  H°(0,1) 

Jo 

with  b  €  H°((0,1)  x  (0,1))  and  v(x)  =  J*b(x,y)dy.  The  associated  sesquilinear 

form  is  given  by 

i 

°2(q)(<P,<W  =  <(vl  -  G)D<p,Di/j>h 

which  is  readily  seen  to  be  H-semicoercive  but  not  V-coercive.  The  generator 
(3.8)  for  the  semigroup  has  domain 

dom  /!(q)  =  {(<p»  e  V  *  H  |  (1  e  V,  EIDVO  =  0, 

[D(EID2<p)  +  (vl  -  G)D0](1)  =  0). 


I 


Time  hysteresis  damping:  These  models,  referred  to  as  Boltzmann  damping 
models,  have  been  widely  studied  in  recent  years  (see  e.g.  [F],  [BF],  [HW]  and 
the  references  therein)  in  connection  with  flexible  structures.  In  an  Euler- 
Bernoulli  beam  as  formulated  above,  the  damping  term  B(q)u(t)  can  be 
replaced  by  a  term  that  has  the  form 


-r 


g(s)D2u(t+s,  )ds 


where  q2  =  g(s)  =  a e^'/^-T,  a,8  >  0.  This  means  that  equation  (3.1)  becomes  a 
functional-partial  differential  equation  or  partial  differential  equation  with 
delay  or  hereditary  term.  Thus  the  above  theoretical  framework  is  not 
directly  applicable.  However,  as  is  shown  in  [BFW],  the  framework  developed 
in  this  section  can  be  appropriately  modified  and  extended  to  give  a  succinct 
theoretical  treatment  of  approximation  methods  involving  estimate 
convergence  and  continuous  dependence  on  observations  for  these  models  also. 


We  have  used  methods  based  on  the  theoretical  ideas  in  this  section  to 
successfully  estimate  damping  in  a  number  of  flexible  structure  experiments. 
To  date  we  have  studied  viscous,  Kelvin-Voigt,  and  Boltzmann  damping  in 
vibrations  of  composite  beams  (BWIC),  [BFW], 


Finally,  wc  note  that  the  ideas  presented  in  this  paper  can  be  extended  to 


provide  a  framework  for  the  treatment  of  nonlinear  distributed  parameter 
systems.  The  cocrcivity  conditions  are  replaced  by  monotonicity  assumptions 
and,  of  course,  nonlinear  semigroups  (evolution  systems)  and  a  nonlinear 
Trotter-Kato  approximation  theorem  play  fundamental  roles  in  development  of 
this  theory.  Details  can  be  found  in  a  forthcoming  manuscript  [BRR]. 
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